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ABSTRACT 


Typical  acoustic  propagation  in  shallow  water  environments  is  dominated  by 
bottom-interacting  paths.  The  effects  of  rough  bottom  interfaces  and  sediment  volume 
fluctuations  are  investigated  using  model  simulations  A  numerical  study  of  low-frequency 
('200  Hz)  broadband  pulse  propagation  is  presented  and  several  characterizations  of 
bottom  factors  are  examined.  In  particular,  the  variations  of  the  interface  rms  roughness 
and  a  volume  fluctuation  strength  constant  on  the  time  resolution  of  the  broadband  pulse 
are  investigated.  It  is  shown  that  interface  roughness  is  the  dominant  factor  in  the 
degradation  of  pulse  resolution. 
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1  UMPE  Model  Parameters 


I.  INTRODUCTION 


Over  the  past  few  years,  the  interest  in  littoral,  shallow  water  environments  has 
increased  significantly.  The  success  of  many  naval  systems  in  these  environments  depends 
on  our  ability  to  predict  and  interpret  shallow  water  acoustic  propagation  The  sound 
speed  profile  in  shallow  water  typically  consists  of  an  upper  mixed  layer  below  which  the 
sound  speed  decreases  to  the  bottom.  Typical  shallow-water  environments  are  found  on 
the  continental  shelf  with  approximate  depths  of  200  m.  The  wave  propagation  is 
downward  refracting,  that  is,  bottom-interacting  paths  dominate  in  long-range  propagation 
as  depicted  in  Figures  1  and  Figure  2.  Etter  (1991)  summarized  the  effect  of  the  sea  floor 
as: 

•  forward  scattering  and  reflection  loss  (but  complicated  by  refraction  in  the 
bottom), 

•  interference  and  frequency  effects; 

•  attenuation  by  sediments; 

•  noise  generation  at  lower  frequencies  due  to  seismic  activity; 

•  backscattering  and  bottom  reverberation 

The  purpose  of  this  thesis  is  to  quantify  the  influence  of  bottom  variability  on  the 
resolution  of  pulse  transmissions  in  simple,  shallow  water  environments. 

Rough  bottom  interfaces  can  be  modeled  by  assuming  some  realistic  form  of  the 
spectral  shape.  The  forward  scatter  of  the  acoustic  field  can  then  be  treated 
deterministically.  Fluctuations  of  the  sediment  properties,  such  as  porosity,  density, 
permeability,  and  shear  strength,  determine  the  sound  scattering  from  the  bottom  volume. 
Yamamoto  (1995)  summarized  the  relationship  between  spectral  strength  and  sediment 
type.  We  used  this  information  to  generate  random  two-dimensional  wave-number 
spectra  of  the  velocity  variability  in  order  to  simulate  a  realistic  ocean  sediment  volume. 
These  techniques  are  described  in  more  detail  in  Chapter  II. 
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Typical  sound  speed  profile  in  shallow  water 


Figure  1.  Typical  sound  speed  profile  in  shallow  water. 


Ray  trace  in  shallow  water 


Figure  2.  Ray  trace  in  shallow  water. 


Wave  propagation  in  an  inhomogeneous  medium  is  defined  by  the  wave  equation. 
The  solution  can  be  obtained  by  several  approximate  methods,  for  example,  Wentzel- 
Kramers-Brillouin  (WKB)  approximation,  ray  approximation  or  parabolic  equation 
approximation.  The  last  method  was  used  throughout  this  thesis  study.  Over  the  past  few 
decades,  ocean  acoustic  models  have  been  developed  and  improved  to  analyze  measured 
data  and  simulate  sound  propagation  in  real  ocean  environments.  Modeling  provides  an 
efficient  means  by  which  to  parametrically  investigate  the  projector  (source)  radiated 
sound  pattern  under  the  influence  of  environmental  conditions.  A  perfect  model  is  one  that 
perfectly  represents  reality.  In  practice,  such  a  perfect  model  would  be  as  complex  as  the 
problem  it  is  attempting  to  represent  So,  there  are  some  assumptions  built  into  every 
computer  model  This  thesis  uses  a  version  of  the  PE  model  developed  at  the  University  of 
Miami  known  as  UMPE  (Smith  and  Tappert,  1994). 

In  order  to  investigate  the  effects  of  bottom  interface  roughness  and  volume 
fluctuations,  we  simulate  broadband  pulse  propagation  for  four  different  scenarios  as 
follows: 

•  flat  bottom  with  homogeneous  sediment; 

•  rough  bottom  interface  with  homogeneous  sediment; 

•  flat  bottom  interface  with  realistic  sediment  sound  speed  fluctuations;  and 

•  rough  bottom  interface  with  realistic  sediment  sound  speed  fluctuations. 


3 


n.  BACKGROUND 


A.  THE  ACOUSTIC  MODEL 


Because  propagation  ranges  are  typically  much  larger  than  ocean  depths,  it  is 
appropriate  to  use  a  cylindrical  coordinate  system  for  long  range  propagation.  Customarily 
we  define  z  as  depth,  r  as  range  outward  from  the  origin,  and  <f>  as  the  azimuthal  angle 
{0—>2 n)  defining  direction  along  a  radial  in  the  horizontal  plane  as  depicted  in  Figure  3.  In 
general,  full  3-D  modeling  of  the  propagation  is  seldom  done  due  to  mathematical  and 
computational  complexity.  It  is  often  a  good  approximation  to  the  real  ocean  acoustic 
environment  to  consider  it  azimuthallv  invariant,  i.e  ,  to  assume  that  environmental 
parameters  such  as  sound-speed  profile,  water  depth,  and  bottom  composition  are 
independent  of  azimuthal  angle  to  the  source.  With  this  assumption,  we  concentrate  the 
analysis  in  depth  versus  range. 

As  mentioned  in  the  introduction,  the  PE  method  is  used  in  this  thesis  to  model 
acoustic  propagation.  The  parabolic  equation  method  was  introduced  into  underwater 
acoustics  in  the  early  1970s  by  Hardin  and  Tappert  (1973)  who  devised  an  efficient 
numerical  solution  scheme  based  on  fast  Fourier  transforms.  We  will  derive  a  PE  solution 
in  a  2-D  constant-density  fluid  environment  where  we  have  assumed  azimuthal  symmetry 
and  hence  no  dependence  on  the  (f)  coordinate. 

We  begin  with  the  wave  equation  for  a  constant-density  medium  in  cylindrical 
coordinates  ( r,z )  independent  of  <j> 


,  1  d'-p 

v2/>— y- f  =  o. 

c  ol 

We  now  assume  we  have  a  steady  state,  time-harmonic  acoustic  field 

p(r,z,tnt)  =  p(r,z)e~,m  . 


(1) 

(2) 
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Top  View 


Substituting  this  into  Eq.  (1)  leads  to  the  Helmholtz  equation, 

^+~^  +  ^  +  k2°n2p  =  °  ’  (3) 

or'  r  a4  az 

oo 

where  p(r,z)  is  the  acoustic  pressure,  ku  =—  is  the  reference  wave  number, 

cn 


Mr,z)  =  7— ~  r  is  the  index  of  refraction,  cn  is  the  reference  sound  speed,  and  c(r,z)  is  the 
c\r^z) 

acoustic  sound  speed.  It  is  within  c(r,z)  that  all  features  of  the  environment  are 
represented 

Tappert  (1977)  assumed  the  solution  of  Eq.  (3)  in  terms  of  an  outgoing  Hankel 
function  with  an  envelope  function  y/(r,z). 

p(r,z)  =  (p{r,z)H'!,(kjj  ,  (4) 


where  ip(r,z)  is  assumed  to  be  slowly  varying  in  range.  In  the  far-fteld  ( k0r  » 1),  the 
Hankel  function  can  be  approximated  by 


(5) 


thus  an  alternative  relationship  between  the  acoustic  pressure,  p(r,z j,  and  the  slowly 
modulating  envelope  function,  ip(r,z ),  is 


U 


P{r,z)  =  P0yj-yip(r,z)eik°r 


(6) 


This  is  the  standard  definition  of  the  so-called  "PE  field  function"  ip  scaled  such  that  at 
r  -  Ra ,  |^|  =  1  and  |/?|  =  P0 .  Substituting  (6)  into  (3)  yields 


P2y/  „  ,  Pij/  d~\p  .,/  ,  ,\ 

— —  +  2 ik..  — —  +  — —  +  k'.br  -  lit//  =  0 
3-  "  dr  dz-  "v  1 


(7) 


Now  we  introduce  the  crucial  paraxial  approximation  based  on  the  slowly  modulating 
character  of  ip, 


d2\p  dy/ 

«  2lk°  If 


(8) 
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(9) 


We  then  obtain  the  following  wave  equation. 


2 ik„  —  +  C—^r  +  A,;  («:  -  l)  y  =  0  , 


"  dr 


dz~ 


which  is  the  standard  parabolic  equation.  This  can  be  rewritten  as 


dy  ik0  f  2  ,  1  o 


&  \ 


dr 


2  V 


n~ 


1  + 


K  8z 


V  , 


(10) 


which  is  the  form  first  introduced  to  the  underwater  acoustics  community  by  Tappert 
(1974) 

Equation  (10)  is  a  one-way  wave  equation  which  can  be  solved  by  a  range¬ 
marching  numerical  technique.  The  UMPE  model  applied  in  this  thesis  uses  the  split-step 
Fourier  (SSF)  method  (Flardin  and  Tappert,  1973)  which  solves  the  parabolic  equation  by 
Fourier  transform  techniques.  Define  the  complex  Fourier  transform  pair  as 


y(r.z)  =  J  y{i\k)e‘k::dk:  =  F  '|V(/\A\)]  , 

(11) 

y(i\k:)  =  —[  !  y{i\:)e  lk-zdz  =  k[y{i\z)]  , 

(12) 

where  kz  is  the  vertical  wave  number.  We  may  rewrite  Eq.  (10)  as 

fjr=*.(r»+u,y . 

(13) 

where 

T  1 

Ikldz 2  ’ 

(14) 

represents  the  kinetic  energy  operator  and 

C! 

*§ 

II 

to 

1 

(15) 

represents  the  potential  energy  operator.  This  representation  of  the  operators  as  kinetic 
and  potential  energy  operators  is  especially  insightful  when  one  wishes  to  form  the  ray 
equations  which  have  Hamiltonian  form  (Smith,  1992). 

In  the  PE/SSF  method,  solutions  of  y  are  marched  out  in  range  by 


8 


</(r  +  Ar)  =  e~ik^e-ik^>y{r)  .  .  (16) 

The  operator  Uop  is  simply  a  multiplication  operator  in  z-space  but  Top  contains  cross 
terms  which  couple  different  depths.  However,  Top  is  a  multiplication  operator  in  vertical 
wave  number  space.  It  is  therefore  desirable  to  separate  the  application  of  each  operator, 
one  in  z-space  and  one  in  £-space  (here,  k  is  the  vertical  wave  number).  In  other  words, 
the  PE  field  function  y/  is  specified  at  some  range  r  initially  in  the  z-domain  A 
transformation  is  made  to  the  ^-domain  followed  by  a  multiplication  of  the  £-space 

operator  e  where  for  is  the  corresponding  kinetic  energy  operator  in  the  wave 

number  domain  The  result  is  then  transformed  again  to  the  z-domain  and  is  followed  by  a 
multiplication  of  the  z-space  operator  e  The  final  result  is  the  field  function  at 

r  At\  that  is, 

y{r  +  Ar,z)  =  x  F~'  x  [i^(r,z))|  .  (17) 

This  is  the  general  algorithm  of  the  PE/SSF  method. 

Previously,  we  have  assumed  the  operators  took  the  forms  defined  by  Eq.  (14)  and 
( 1 5)  These  forms  are  commonly  referred  to  as  the  “standard  PE”  and  are  only  one  set  of  a 
number  of  various  operator  forms.  In  our  simulation,  we  will  use  a  higher  order  form 
which  was  introduced  by  Thomson  and  Chapman  (1983)  referred  to  as  the  “wide-angle” 
PE  approximation  (WAPE)  with  operators  given  by 

UW=|- 

-l]  (19) 

We  then  need  to  define  the  initial  conditions  for  the  PE  field  function,  y/(r-0,z).  It 
is  desirable  to  begin  the  calculation  by  specifying  the  source  in  the  ^-domain.  Thomson 
and  Bohun  (1988)  have  shown  that  a  wide  angle  source  can  be  defined  as 


(18) 
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The  calculation  of  the  time  domain  arrival  structure  at  a  given  range  is  an  integral 
part  of  the  analysis  in  this  thesis.  The  pressure  field  in  the  time  domain  is  represented  by 

p{r,z,t)  =  F-\pa(r,z)]  =  YjaP<Xr,z)e,<ot  ,  (21) 

where  pm{r,z)  is  a  single  frequency  component.  The  arrival  time  structure  at  some  fixed 

range  r=R  can  be  computed  using  the  complex  field  p^R^z)  for  many  frequencies  and 

then  Fourier  transforming  to  obtain  p(R,z,t),  the  set  of  complex  pressure  values  in  the 
time  domain  To  avoid  large  transform  sizes,  the  UMPE  model  basebands  the  solutions 
and  shifts  the  center  frequency  bin  to  zero.  To  reduce  the  effect  of  side-lobes,  we  apply  a 
Hanning  filter  over  the  bandwidth  before  transforming  to  the  time  domain.  Additionally, 

R 

ia> — 

we  neglect  the  phase  factor  e  °  -e  c°  by  simply  performing  the  Fourier  analysis  on  the 
field  function  y/a(R,z).  Therefore,  the  time  domain  is  heterodyned  around  the  value 
I,,  -R/Crj  and  arrival  times  are  computed  as  values  of  “reduced  time”  or  (t-t0). 

C.  ROUGH  BOTTOM  INTERFACE 


The  ocean  bottom  roughness  is  assumed  to  be  characterized  by  a  2-D  power  law 
spectrum  at  high  wave  numbers  which  can  be  written  as 


where  Lajrr  is  the  correlation  length  of  the  roughness  and  /?  is  the  spectral  exponent.  (Note 
that  our  use  of  wave  number  is  now  with  respect  to  range  scales  of  interface  roughness.) 
To  give  the  full  spectrum  a  realistically  smooth  structure,  we  assume  the  spectral  form 
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(23) 
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The  roughness  realization  is  then  simply  rescaled  by  the  desired  value  of  o  to  provide  the 
appropriate  rms  value. 

D.  SEDIMENT  COMPOSITION 


Yamamoto  (1995)  summarizes  the  transformation  between  three-dimensional 
wave  number  spectra  and  the  velocity  variabilities  which  determine  the  ocean  sediment 
properties  This  transformation  can  be  simplified  for  our  2-D  case.  The  2-D  compressional 
wave  velocity  structure  in  the  sediment  may  be  written  as 

c(x,z)  =  c^(l  +8c(x,zj)  (30) 

where  c(x,z)  is  the  velocity  at  a  point  (x,z),  c0  is  the  average  velocity  at  the  seafloor  and 
5c(x,z)  is  the  relative  variability.  The  relative  variability  is  modeled  as  a  zero-mean 
quasistationary  random  function,  i.e., 

(4(x,r))  =  0  (31) 

(Sc(x,z)Sc(x',y'))  =  [  J  /  S{kx,k_)xcos{kx{x-x’)  +  kXz-z')}dkxdk_  (32) 
The  2-D  power  spectrum  of  the  relative  velocity  variability  is  denoted  by  S(kx,kv) ■  By- 
definition,  the  variance  (S  represents  the  total  energy  of  the  power  spectrum. 


fe:}=r <33) 

It  is  a  plausible  assumption  that  the  vertical  power  spectrum  may  be  given  by  a  power-law 
form 


S,(k)=Bk;r  (34) 

where  /?  is  the  spectral  exponent  and  B  is  the  spectral  strength  constant.  The  2-D  power 
spectrum  given  by  Yamamoto  is 


$M.-) 


where  a  is  the  aspect  ratio  and  /?  2  must  be  an  integer 


(35) 
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To  produce  a  random  realization  of  sediment  volume  fluctuations,  we  assume  the 
complex  wave  number  spectrum  has  a  random  amplitude  and  phase,  Ae'^,  for  each  wave 
number  component.  The  random  amplitude  of  each  component  of  the  power  spectrum  can 
then  be  obtained  by  A2  =-ln(r)  where  r  is  an  independent,  uniformly  distributed  random 
variable  in  the  interval  (0,1).  The  random  phase  <j)  is  another  independent,  uniformly 
distributed  random  variable  in  the  interval  (0,27t).  To  avoid  infinite  magnitudes  at  low 
wave  numbers,  the  spectrum  is  high-pass  filtered.  The  filtered  spectrum  is  then  Fourier 
transformed  to  obtain  a  realization  of  the  sound  speed  variability  The  entire  field  of 
sediment  fluctuations  is  first  computed  on  a  fine  grid  and  then  interpolated  to  match  the 
computational  grid  of  the  PE/SSF  algorithm  The  entire  scheme  can  be  described  as 
follows: 


Step  1.  Construct  the  2-D  wave  number  field  (kr  (n),  k,  (w)) ,  according  to 


M„  = 


in 


kr{n)-{n-\)Akr  n  =  1,2,..., N  , 


(36) 


A .  2n 

A  [z  -z  )  ’ 

V^max  Abot ) 

k2{m )  =  (m~  OA^  m  =  1,2, . . . ,  M  -  Mbot  , 

where 

rmax  =  maximum  propagation  range, 

N  =  FFT  size  in  range  direction, 
jm;iX  =  maximum  compution  depth, 

:hijl  ~  average  bottom  depth 
M  -  FFT  size  in  depth  direction, 

MhiH  -  index  of  average  bottom  interface. 


(37) 


Step  2.  Generate  two  random  numbers  series  for  the  random  phase  and 
amplitude, 

I - 7 - C 

(38) 


A(n,m)  =  }/-ln(anm)  , 


13 


<t>{ri,m)  =  2 jtbnm  ,  (39) 

where 

a,b  =  uniform  random  number  distribution  0  -»  1 . 

•  Step  3.  Compute  the  random  phase  and  amplitude  for  every  mesh  point  in  the 
2-D  realized  wave  number  field,  S'(kr(n)k.(m)} ,  i.e., 

S'(k, .(//), kirn))  -  S(kt  {n)X{m))A{njv)emnm]  .  (40) 

•  Step  4.  Add  high-pass  filter  to  normalize  the  exponent  amplitude  at  low  wave 
number 

•  Step  5.  Compute  2-D  FFT  to  yield  relative  velocity  variability, 

Sc(r,z)  =  2 DFFT\S'{kr{n)k2{m))\ . 

The  velocity  structure  of  the  sediment  realization  is  then  given  by: 

c(r,  z )  =  cu  +  C(t  x  R e(S£  (r,  zfj  (41) 

E.  BEAMFORMING 

In  underwater  acoustic  propagation,  much  of  the  information  about  the 
environment  and  the  relative  source/receiver  locations  is  contained  in  the  vertical  structure 
ot  the  acoustic  field.  This  is  especially  true  in  shallow  water  environments  where  the 
vertical  structure  can  be  defined  in  terms  of  a  relatively  small  number  of  modes  of  the 
acoustic  waveguide.  Because  these  modes  are  generally  dispersive,  the  ability  to  resolve 
the  vertical  structure  of  the  time  arrival  at  the  receiver  can  lead  to  estimations  of  the 
corresponding  source  location. 

Compared  to  the  directivity  of  single  sensors,  a  beamformer  permits  one  to  listen 
preferentially  to  wavefronts  propagating  from  one  direction  over  another.  In  this  sense,  a 
beamformer  implements  a  spatial  filter  (see,  for  example,  Defotta  et  al. ,  1988,  Ziomek, 

1 995).  Consider  a  line  array  of  equally  spaced  elements  positioned  along  the  y-axis  in  an 
infinite,  homogeneous  medium  with  an  undetermined  number  of  remote  sources.  Let  the 
output  of  an  element  located  at  the  origin  of  coordinates  due  to  the  Ith  source  be  s/(t). 
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Under  the  assumption  of  plane- wave  propagation  in  the  far-field,  a  source  from  direction 
6,  produces  the  following  sensor  outputs: 


(  nd  sin(&P 

en(t)  =  se  *  + - 7^  .  (4: 

v  c  ) 

where  n  is  the  array  element  index,  d  is  the  spacing,  0t  is  the  angle  of  arrival  from  the  Ith 

source,  and  c  is  the  speed  of  propagation.  Beamforming  is  accomplished  by  applying 
weights  and  time  delays  to  the  individual  element  signals,  and  then  coherently  summing 
them  The  beamformer  output  when  viewed  in  the  direction  0m  is  then 

a  '  f  //j(sin(#f )  -  sin(#m))l 


hmU)  =  Z  w»sf  t 


where  the  w'„  are  weights  applied  to  each  element’s  signal. 

The  directional  response  characteristics  of  the  beamformer  can  be  derived  by 
considering  a  complex,  plane-wave  cw  signal  propagating  across  the  array, 

s,{t)  =  e^  .  (44) 

The  source  produces  the  following  sensor  outputs  due  to  the  arrival  angle 


«„(<)  = 


'.k{d  sm(  8^ 


where  k(  =  =  —  is  the  wave  number  of  the  Ith  source.  The  beamformer  output  is  then 


bjj )  =  =  w(ef)e‘a 


where 

W(6f )  =  w ne'nk>‘"™<"'  ’  s"’^» ,! ,  (47) 

n=  0 

represents  the  spatial  Fourier  transform  of  the  array  weights. 

A  two-dimensional  FFT  (spatial  and  temporal)  can  be  used  to  implement  the 
beamforming  operation  for  an  equally  spaced  line  array.  For  example,  consider  a  line  array 
having  N-8  equally  spaced  elements  as  depicted  in  Figure  4.  The  signals  from  each 
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element  are  sampled  simultaneously  at  rate  fs  to  obtain  a  multichannel  time  series  x„.„.  M- 
point  FFTs  are  taken  along  each  channel  to  obtain  the  temporal  frequency  response,  Xn(/J, 
where  y  is  the  frequency  bin  number,  i.e.. 


\i  i 

A'»(./)  =  X  >TT,,,T  ,  o  <  y  <  M  -  1  .  (48) 

nr  0 

The  Xn(j)  are  represented  in  the  figure  as  N  horizontal  rows  of  M  FFT  bins. 

Next,  the  N  cross-channel  FFT  is  performed  with  appropriate  weighting  along  the 
array,  i.e., 

=  ,0Si<Af-l  (49) 

n=ft 

The  result  of  the  N  cross-channel  FFT  series  to  spatial  frequency  (or  wave  number) 
provides  the  phase  difference  information  in  depth  for  specific  temporal  frequencies.  The 

relationships  between  the  various  indices  and  the  variables  used  in  the  beamformer  output 
are 


/,= 

K  - 


k  2k 
N  d  ’ 


(50) 


(51) 


=  sin“ 


c  kt 

{2*  ft) 


(52) 


wherey  and  k  are  the  indices  in  the  frequency  and  wave  number  domains,  respectively.  In 
practice,  6m  is  determined  by  interpolating  between  values  of  ki .  After  the  data  has  been 
interpolated  to  arrival  angle/frequency  space,  an  inverse  FFT  is  taken  along  each  angle  bin 
to  provide  arrival  angle  as  a  function  of  travel  time. 
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Figure  4.  Beamforming  example. 
From  [DeFotta,  1988] 
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IE.  MODEL  SIMULATIONS 


As  was  described  in  the  introduction,  the  principal  characteristic  of  most  shallow- 
water  propagation  is  dominated  by  bottom-interacting  paths  Sound  transmission  in 
shallow  water  has  an  optimum  frequency  of  propagation  at  longer  ranges  as  shown  by 
Jensen  and  Kuperman  (1983).  Optimum  frequency  is  a  general  feature  of  ducted 
propagation  in  the  ocean.  It  occurs  as  a  result  of  competing  propagation  and  attenuation 
mechanisms  at  high  and  low  frequencies.  In  the  high-frequency  regime  there  is  increasing 
volume  and  scattering  loss  with  increasing  frequency.  At  lower  frequencies  the  situation  is 
more  complicated.  With  increasing  wavelength  the  efficiency  of  the  waveguide  to  confine 
sound  decreases  (the  cutoff  phenomenon).  Hence  propagation  and  attenuation 
mechanisms  outside  the  waveguide  (in  the  seabed)  become  important.  In  fact,  the 
increased  penetration  of  sound  into  a  lossy  seabed  with  decreasing  frequency  causes  the 
overall  attenuation  of  waterborne  sound  to  increase  with  decreasing  frequency.  Thus  we 
get  high  attenuation  at  both  high  and  low  frequencies,  while  intermediate  frequencies  have 
the  lowest  attenuation.  As  Jensen’s  study  shows,  the  optimum  frequency  is  in  the  range 
200-800  Hz  for  a  water  depth  of  200  m  and  we  will  use  200  Hz  as  our  center  frequency 

We  assume  that  the  ocean  bottom  sediment  is  a  fluid-like  half-space  which  means 
that  it  supports  only  one  type  of  sound  wave,  i.e.,  compressional  waves.  The  information 
required  for  model  simulations  should  include  the  following  depth-dependent  material 
properties:  the  sound  speed  profile  in  the  water  column  and  compressional  wave  speed, 
attenuation  and  the  density  in  the  seabed.  For  all  cases,  we  will  assume  the  water  column 
is  homogeneous  with  a  constant  sound  speed. 

This  chapter  presents  the  model  simulation  results  of  four  test  cases.  Case  1 
examines  the  propagation  of  a  broadband,  impulsive  source  over  a  flat,  homogeneous 
bottom.  Case  2  examines  the  propagation  in  a  similar  environment  with  a  rough  bottom 
interface.  We  then  examine  which  parameters,  i.e.,  rms  roughness,  correlation  length  or 
ratio  of  them,  dominates  the  beam  resolution.  Case  3  is  also  similar  to  Case  1  but  with 


18 


realistic  bottom  volume  fluctuations  The  factors  which  most  strongly  influence  resolution 
in  Case  2  and  Case  3  will  be  combined  into  Case  4 

The  basic  range-independent  parameters  used  in  the  UMPE  model  are  listed  in 
Table  1. 


Model  Type 

Wide  Angle  PE 

Source  Function  Type 

Wide  Angle  Source 

Reference  Sound  Speed 

1500.00  m/s 

Bottom  Depth 

200.0  m 

Maximum  Computational  Depth 

400.0  m 

Mesh  Points  For  Positive  Depths 

512 

Maximum  Range 

10.0  km 

Step  In  Range 

2.0  m 

Mesh  Points  In  Range 

5000 

Source  Depth 

100.0  m 

Center  Frequency 

200.0  Hz 

Frequency  Bandwidth 

100.0  Hz 

Number  Of  Frequencies 

256 

Sound  Speed  In  Water 

1500  m/s 

Sound  Speed  In  Sediment 

1600  m/s 

Density  In  Sediment 

1 500  kg/m3 

Attenuation  In  Sediment 

0.2  dB/m/kHz 

Table  1.  UMPE  Model  Parameters 


A.  CASE  1:  RANGE  INDEPENDENT 

We  begin  with  the  simple  case  consisting  of  a  constant  speed  water  column  and 
fluid-like  homogeneous  sediment  separated  by  a  flat  interface  The  corresponding 
transmission  loss  field  for  the  center  frequency  is  depicted  in  Figure  5.  We  note  the 
symmetric  constructive  interference  due  to  the  finite  number  of  trapped,  propagating 

normal  modes.  By  Snell's  law,  there  is  a  critical  angle  at  arc  cos  *  ~  ^  =20.4°  (relative  to 

1600 

horizontal)  below  which  no  energy  is  transmitted  into  the  bottom.  From  the  UMPE  output 
field  with  beamforming  at  10,000  m,  we  obtain  the  angle-time  plot  shown  in  Figure  6.  As 
expected,  this  energy  arrives  primarily  near  horizontal  with  an  aperture  of  approximately 
±20°.  The  steeper-angle  energy  ^20°  is  rapidly  lost  by  transmission  into  the  seabed. 
Furthermore,  the  lowest  modes  arrive  first,  before  the  higher  modes,  as  would  be  expected 
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Case  1  (Flat  bottom) 


(tu)  ijictea 


Figure  5.  Range-depth  transmission  loss  field  for  flat  bottom  interface. 
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Figure  6.  Arrival  angle-time  plot  for  flat  bottom  interface. 


from  an  analysis  of  the  modal  dispersion  in  this  environment  A  standard  calculation  of  the 
number  of  trapped  modes  was  found  to  be  consistent  with  the  number  of  arriving  modal 
peaks  when  the  range  of  the  transmission  loss  palette  was  extended  to  expose  the  weakest 
modal  arrivals. 

B.  CASE  2:  ROUGH  BOTTOM  INTERFACE 

A  rough  seafloor  will  cause  attenuation  of  the  mean  acoustic  field  propagating  in 
the  ocean  waveguide  due  to  incoherent  scattering  thereby  decreasing  the  mean  signal  to 
noise  ratio.  For  individual  pulse  transmissions,  the  instantaneous  signal  becomes  smeared 
in  time  and  space  due  to  mode  coupling.  A  rough  bottom  interface  is  added  to  the 
environment  of  Case  1  in  order  to  investigate  the  effect  of  scattering  on  the  instantaneous 
resolution  of  the  pulse  in  the  water  column. 

As  described  in  Chapter  II,  there  are  two  parameters  that  control  the  roughness, 
i.e,,  the  rms  roughness  and  the  correlation  length  In  order  to  determine  the  dominant 
factor,  we  simulated  two  sets  of  comparisons.  In  one,  the  correlation  length  is  fixed  at 
1 000  m  and  the  rms  roughness  has  values  of  2,  5,  and  1 0  m.  In  the  other,  the  rms 
roughness  is  fixed  at  2  m  with  correlation  lengths  of  40,  200,  and  400  m.  Figure  7-12 
show  the  beamformed  angle-time  plots  at  10,000  m.  We  can  see  that  increasing  the  rms 
roughness  and  decreasing  the  correlation  length  will  reduce  the  resolution  but  it  is 
difficult  to  see  which  factor  dominates.  However,  the  distortion  changes  look  similar  and 
suggest  that  the  ratio  of  rms  roughness  and  correlation  length  may  be  the  important 
parameter.  We  simulated  several  realizations  with  fixed  ratio  of  the  rms  roughness  and  the 
correlation  length  to  confirm  our  hypothesis,  specifically,  (rms  roughness):  (correlation)  = 
1:500,  1:100  and  1:20.  The  beamformed  results  are  displayed  in  Figure  13-15.  Combined 
with  the  previous  simulations,  we  find  that  higher  ratios  retain  higher  resolution  and  lower 
ratios  degrade  the  resolution  more  severely.  Figure  16-18  show  single  frequency  (200  Hz) 
transmission  loss  fields  in  depth-range  slices  for  various  rough  interface  realizations. 
Compared  to  the  flat  bottom  case,  we  see  that  there  is  more  energy  penetration  into  the 
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Figure  7.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
2m  and  correlation  length  =  1,000m. 
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Case2  (rms=5m,  correlation  length= 1,000m) 


o 

*6 


o 

d 

CD 


o 

o 

o 

o 

o 

d 

d 

o 

o 

d 

CO 

CO 

CO 

CO 

(Sap)  ajStre  p^Aixiy 


Figure  8.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
5m  and  correlation  length  =  1,000m. 


TL  (dB) 


(Sap)  ajStre  [BAUjy 


Figure  9.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
10m  and  correlation  length  =  1,000m. 
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Figure  10.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
2m  and  correlation  length  =  40m. 
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Figure  11.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
2m  and  correlation  length  =  200m. 
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Figure  12.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
2m  and  correlation  length  =  400m. 
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Figure  13.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
5m  and  correlation  length  =  2,500m. 
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Figure  14.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
2m  and  correlation  length  =  200m. 
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Figure  15.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
5m  and  correlation  length  =  100m. 
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Figure  16.  Range-depth  transmission  loss  Held  for  rough  bottom  interface  with  rms 
roughness  -  5m  and  correlation  length  =  2,500m. 
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Case  2  (rms=2m,  corr.=200m) 
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Figure  17.  Range-depth  transmission  loss  field  for  rough  bottom  interface  with  rms 
roughness  =  2m  and  correlation  length  =  200m. 
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(UJ)  md0Q 


E 


CD 

O) 

c 

CO 

QC 


Figure  18.  Range-depth  transmission  loss  field  for  rough  bottom  interface  with  rms 
roughness  =  5m  and  correlation  length  =  100m. 
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bottom  at  longer  ranges,  and  the  stable  interference  structure  has  broken  down  into  a 
complex  pattern. 

C.  CASE  3:  SEDIMENT  VOLUME  FLUCTUATIONS 

The  scattering  due  to  the  velocity  variability  within  the  sediment  is  considered  to 
be  the  major  scattering  mechanism  in  many  geological  fields  in  shallow  water  based  on 
extensive  measurements  of  acoustic  backscattering  and  sediment  cores  (Jackson  and 
Briggs  1992,  Yamamoto  1995).  The  geology  of  shallow-water  seabeds  is  complicated 
because  of  the  physical,  chemical,  and  biological  activities  which  took  place  though  the 
long  geological  time  spans.  The  measurements  of  the  velocity  and  density  variabilities  in 
the  ocean  are  beyond  this  thesis  subject.  For  our  simulations,  we  shall  use  the  sediment 
model  described  in  Chapter  II  with  some  typical  values  extracted  from  Yamamoto's 
work. 

In  this  case,  we  use  a  flat  bottom  interface  and  add  the  sediment  volume 
perturbation.  According  to  Yamamoto's  model,  the  bottom  sediment  perturbation  is 
controlled  by  three  parameters:  spectral  exponent  P,  aspect  ratio  a,  and  the  spectral 
strength  constant  B.  In  his  conclusions,  the  spectral  exponent  P  and  aspect  ratio  a  are 
fairly  constant  regardless  of  the  sediment  type  but  the  spectral  strength  constant  B  varies 
depending  on  the  sediment  type  and  location.  We  will  fix  a  and  P  and  vary  the  values  of 
B  by  choosing  three  typical  spectral  strength  constants,  e.g.,  2x10 3,  2x10  4  and  2x10 5 
(units  of  B  are  aw1'*5).  The  sediment  field  for  strength  constant  B  =  2xl0‘3  wl'p  is  depicted 
in  Figure  19.  The  rms  value  of  the  sound  speed  fluctuations  in  the  sediment  for  this  value 
of  B  is  22.4  m/s.  The  beamformed  results  for  the  three  values  of  B  are  shown  in  Figure 
20-22.  The  pulse  shape  does  appear  to  degrade  slightly  as  the  spectral  strength  constant  is 
increased.  Overall,  however,  the  sediment  volume  perturbation  doesn't  seem  to  affect  the 
pulse  resolution  significantly. 
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Sediment  variability  (B=2e-3) 
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Amplitude  of  velocity  variability  (m/s) 
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Figure  20,  Arrival  angle-time  plot  for  sediment  volume  fluctuation  with  strength 
constant  B-2xl0  3ml  p. 
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Figure  22.  Arrival  angle-time  plot  for  sediment  volume  fluctuation  with  strength 
constant  B=2xl0':,m1'p. 
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D.  CASE  4:  ROUGH  BOTTOM  INTERFACE  WITH  SEDIMENT  VOLUME 

FLUCTUATIONS 

We  have  seen  the  effect  of  rough  bottom  interfaces  and  sediment  volume 
fluctuations  separately.  We  now  need  to  combine  these  two  factors.  From  Case  2  and 
Case  3  results,  it  appears  that  the  effect  of  a  rough  bottom  interface  dominates  the 
pulseresolution.  The  strongest  sediment  volume  perturbation  strength  constant,  B  =  2x1  O'3 
/n1'15,  was  chosen  to  produce  the  largest  volume  effect.  The  interface  correlation  length 
was  fixed  at  1000  m  and  the  rms  roughness  was  varied  from  2  m  to  10  m.  Comparing 
these  results  to  those  cases  with  no  volume  fluctuations  showed  little  difference  indicating 
that  these  levels  of  interface  roughness  dominate  the  effects  on  pulse  resolution. 

E.  PULSE  RESOLUTION  ANALYSIS 

We  have  shown  that  the  effect  of  the  bottom  roughness  and  sediment  volume 
fluctuation  degrades  low  frequency  pulse  resolution  in  shallow  water  propagation.  In  this 
section,  we  will  quantify  the  effect  of  rough  bottom  and  sediment  volume  fluctuation  on 
pulse  resolution  in  order  to  obtain  a  clearer  relationship.  The  conventional  angle  resolution 
is  defined  as  an  angle  width  of  Me  peak  amplitude.  In  order  to  simplify  the  Case  2  analysis, 
we  fixed  the  correlation  length  at  a  value  of  1000  m  and  varied  the  rms  roughness.  Figure 
23  shows  several  slices  in  reduced  times  of  1,  1.5,  2,  2.5  and  3  seconds  from  an  arrival 
angle-time  plot  for  rms  roughness  of  0  (the  Case  1  results),  2,  5,  10  and  20  m.  Each 
reduced  time  data  point  is  actually  an  average  of  five  points  spread  over  a  time  bin  width 
of  0.2  sec,  i.e.,  the  datum  for  the  reduced  time  of  1  sec  is  an  average  of  data  at  0.90,  0.95, 
1.00,  1.05,  and  1.10  sec.  Polynomial  fits  to  the  data  produce  the  smooth  curves  shown  in 
the  figures.  We  can  see  a  good  angle  resolution  at  lower  rms  roughness  values.  As  the  rms 
roughness  is  increased,  the  resolution  is  found  to  degrade,  particularly  for  the  later  arrival 
times  This  is  expected  since  later  arrivals  correspond  to  higher  modes  which  interact  more 
strongly  with  the  bottom.  As  an  estimation  of  this  effect,  we  note  that  the  angular  pulse 
width  appears  to  double  for  an  rms  value  of  ~10  m  corresponding  to  an  rms:  correlation 
ratio  of  1: 100. 
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In  Figure  24.  slices  are  taken  at  arrival  angles  0,  5,  10,  and  20  deg  In  this  case, 
three  points  over  a  width  of  1  deg  are  averaged  to  produce  each  data  point  shown.  There 
is  an  anomalous  point  at  an  rms  of  2  m  due  to  unresolved  peaks  of  two  modes.  This  point 
was  not  used  in  the  curve  fit.  The  time  resolution  degrades  when  the  rms  roughness 
increases  for  all  arrival  angles  except  0  deg.  This  is  due  to  less  bottom  interactions  of  the 
horizontal  propagation.  Furthermore,  we  notice  that  the  time  width  does  not  degrade  as 
rapidly  as  the  angular  width.  None  of  the  rms  roughness  values  appear  to  double  the  pulse 
width  in  time.  T  he  peak  minimum  transmission  loss  also  increases  with  rms  roughness  as 
shown  in  Figure  25  by  roughly  5  -  10  dB  at  an  rms  value  of  10  m. 

To  perform  a  similar  analysis  of  the  Case  3  environment,  the  spectral  exponent  and 
aspect  ratio  were  fixed.  Obviously,  the  sediment  volume  fluctuations  show  more  stable 
results  with  varying  strength  constant  as  depicted  in  Figure  26-28.  The  angle  and  time 
resolution  change  little  with  increasing  strength  constant.  The  peak  (minimum) 
transmission  loss  also  remains  stable. 

In  Figure  29,  the  peak  minimum  transmission  loss  computed  with  the  combined 
effects  of  a  rough  bottom  interface  and  sediment  volume  fluctuations  shows  very  similar 
results  as  in  Figure  25.  It  is  reasonable  to  conclude  that  these  levels  of  bottom  roughness 
dominate  pulse  resolution  effects  in  simple  shallow  water  propagation. 
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Figure  23.  Angle  width  versus  rms  roughness  for  reduced  time  =  1,  1.5,  2,  2.5  and  3 
second. 
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Figure  24.  Time  width  versus  rms  roughness  for  arrival  angle  =  0,  5, 10, 15  and  20 
degree. 
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Figure  26.  Angle  width  versus  strength  constant  for  reduced  time  =  1, 1.5,  2,  2.5  and 
3  second. 
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Figure  28.  Peak  minimum  transmission  loss  versus  strength  constant  for  reduced 
time  =  1,  1.5,  2,  2.5  and  3  second. 
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Figuie  29.  Peak  minimum  transmission  loss  versus  rms  roughness  for  reduced  time 
=  I,  1.5,  2,  2.5  and  3  second  combine  with  strength  constant  =  2x10  3  m1  p  . 
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IV.  CONCLUSION 


In  this  thesis,  we  have  presented  a  numerical  investigation  of  the  effects  of  bottom 
interface  roughness  and  sediment  volume  fluctuations  on  the  resolution  of  broadband 
pulse  propagation  in  a  simple,  shallow  water  environment.  Predictions  of  the  acoustic 
propagation  were  accomplished  with  a  broadband,  wide-angle  parabolic  equation  model. 
Fluctuations  of  the  bottom  interface  and  volume  were  computed  from  realistic  spectra 
characterizing  each  structure.  The  strength,  or  level  of  variability,  of  each  was  varied 
independently  in  order  to  determine  the  relative  influence  on  pulse  resolution 

To  be  consistent  with  measurement  techniques  employed  in  underwater  acoustics 
experiments,  a  beamforming  algorithm  was  written  to  discriminate  vertical  arrival  angles 
as  a  function  of  arrival  time.  When  no  fluctuations  were  introduced,  the  solution  to  the 
range-independent  problem  was  found  to  agree  with  a  standard  normal  mode  analysis. 
When  fluctuations  were  introduced,  the  symmetrical  arrival  time  pattern  became  distorted, 
indicating  a  degradation  in  the  resolution  of  the  pulse.  Such  degradations  are  expected  to 
limit  the  ability  of  some  operational  naval  systems 

When  interface  roughness  was  introduced,  the  results  suggested  that  the  main 
parameter  influencing  resolution  was  the  ratio  of  rms  roughness  to  correlation  length.  In 
other  words,  a  high  value  of  rms  roughness  does  not  independently  imply  a  degradation  in 
the  pulse  shape  if  the  correlation  length  is  quite  large.  Similarly,  a  relatively  small  rms 
roughness  may  still  significantly  affect  pulse  resolution  if  accompanied  by  a  very  small 
correlation  length  In  contrast,  the  largest  realistic  scales  of  volume  fluctuations  had  little 
impact  on  the  pulse  shape.  A  combination  of  large  volume  fluctuations  and  various  scales 
of  rms  interface  roughness  showed  that  degradations  of  pulse  resolution  are  dominated  by 
the  interface  structure.  While  it  may  be  possible  that  a  different  environment  (e.g.,  strongly 
downward  refracting  water  column  and  upward  refracting  sediment)  would  increase  the 
influence  of  sediment  heterogeneity,  such  cases  would  also  tend  to  increase  interactions 
with  the  bottom  interface  resulting  in  similar  conclusions. 
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An  attempt  was  made  to  quantify  these  results  by  measuring  the  1/e  -  width  of  the 
pulse  at  various  arrival  angles  or  arrival  times  as  a  function  of  the  fluctuation  strength.  It 
was  observed  that  an  rms  interface  roughness  of  10  m  with  a  correlation  length  of  1000  m 
produced  a  doubling  of  the  angular  width  of  the  arrivals  and  increased  the  transmission 
loss  by  as  much  as  10  dB.  The  time  width  was  more  stable,  however.  In  none  of  the  cases 
with  volume  fluctuations  was  the  signal  found  to  degrade  appreciably.  When  the  two 
influences  were  combined,  the  results  were  nearly  identical  to  the  case  with  no  volume 
fluctuations. 

Overall,  these  results  suggest  that  bottom  interactions  in  shallow  water 
propagation  are  dominated  by  interface  scatter.  The  degradation  of  pulse  resolution 
depends  on  a  combination  of  the  rms  interface  roughness  and  an  associated  correlation 
length  scale.  In  order  to  know  a  priori  the  limitations  of  pulse  resolution,  it  is  necessary  to 
obtain  an  accurate  statistical  representation  of  the  sea  floor  topography.  In  contrast,  the 
sub-bottom  volume  fluctuations  do  not  appear  to  significantly  influence  the  forward 
propagation 
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